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We consider an "oriented" chiral condensate produced in the squeezed 
states of the effective field theory with time- and space-dependent pion mass 
parameter. We discuss the general properties of the solution, identifying 
condensate modes and determining the resulting pion distributions. The im- 
plementation of the dynamics in the form of sudden perturbation allows us to 
look for exact solutions. In the region of condensation, the dramatic increase 

in pion production and charge fluctuations are demonstrated. 
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§5 . I. INTRODUCTION 

Pions have always played an important role in heavy ion reactions. With their relatively 
strong coupling to nuclei and very small mass, pions can be produced in large quantities 
and can carry important detectable information about the state of nuclear matter during 
the reaction. Large fluctuations in the ratio of charged to neutral pions have been observed 
in cosmic ray experiments |I[ and have recently been a subject of many discussions. The 
^ major interest in this direction is motivated by an almost a decade old idea @-§] that at 

high enough energies the pion-sigma ground state symmetry breaking condensate can be 
destroyed. In the latter process of sudden cooling there is a chance of the formation of 
disoriented chiral domains, given that the terms responsible for the explicit breaking of 
chiral symmetry are small. The formation of this disoriented chiral condensate (DCC) is 
presumably responsible for the observed discrepancy in the ratio of pion species. Further 
enthusiasm about these ideas is generated by the fact that it may be possible to achieve 
an environment capable of producing disoriented chiral domains in relativistic heavy ion 
collisions. In that case the large fluctuations in pion types can be a signal of a DCC. 

There exist a number of theoretical works investigating the production and development 
of DCCs in heavy ion collisions. The assumption of a random equiprobable orientation 
of the pionic isovector leads [[| to the probability of observing a fraction of neutral pions 
/ = N w o/(N w+ + N w - + N w o) being 
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This result is consistent with the slowly varying classical pionic field as a solution of the 
non- linear a-model 0. It can also be obtained as a limiting large N distribution in the 
coherent single-mode pion production by an isoscalar operator 

The formation and the nature of domain structures is an important question in itself P1,|T0 



In the case of many small domains or in the absence of a DCC, one would expect a Gaussian 
distribution P(f) , following from the central limit theorem. More complicated methods 
of quantum and classical field theory applied to this problem with different assumptions of 



formation, evolution [TT].|l2| and dissociation of DCCs [jlcjpl lead to different results. Pionic 
mode mixing and final state interaction may also greatly change the observed forms of these 
probability distributions [§||[15|. 

Most of the existing theoretical works consider pion production from DCC domains 
formed in the dynamical process of symmetry restoration and later sudden relaxation into a 



non-symmetric vacuum which is some times called quench [|T^] . The word "disoriented" also 



describes the theoretical foundation of the DCC approach which is based on introducing an 
effective disoriented current-type term in the Lagrangian. The advantage of this formalism 
is in the simplicity of the solution, as a current-type interaction can always be solved with a 



coherent state formalism |T7[]. In contrast to that the linear sigma model, a well-established 
effective theory based on chiral symmetry, naturally leads to a quadratic pion field term in 
the Lagrangian. The goal of the present work is to consider an "oriented" chiral condensate 
with a chirally symmetric Lagrangian that to the lowest order is quadratic in the pion field. 
Throughout all this paper we are going to deal with the wave equation for pions in the form 

_Z_V 2 7? + m 2 ff (x,£)7r = 0. (2) 

This is the equation of motion deduced from a quadratic Lagrangian. Equation @ contains 
an effective mass m e g (x, t) which is due to a mean field from non-pionic degrees of freedom. 
This term produces a parametric excitation of the pionic field, and, if m^ < 0, may lead to 
amplification of low momentum modes and the formation of a chiral condensate. The chiral 
condensate in this picture is "oriented"; as opposed to the usual picture of a DCC. The 
scalar mean field (effective mass) introduced here does not break the symmetry properties 
of pion fields in any way. The quantum problems formulated by equations of the form @ 
address a variety of physical issues, and can be encountered in different areas of physics. In 
general, such non-stationary quantum problems cannot be solved exactly. Even studies with 
approximate methods like perturbation theory, impulse approximations, adiabatic expan- 
sions and many others often require sophisticated approaches. We will address those rare 
occasions when exact solutions may be obtained. The significance of the analytical solution 
should not be underestimated. The exact results could exhibit the unperturbative features 
of the solutions, like phase transitions and condensates, as well as point the way to a good 
approximate theory. 

The paper is structured as follows. We start with a short introduction to the linear 
sigma model and show that effectively it leads to Eq. (^) for a pion field. In the next 
section we discuss an exact solution of Eq. (0). First, we show that from the classical 
solution a quantum solution can be built exactly. Then, based on the general form of the 
evolution matrix, we obtain and discuss possible forms of the multiplicity distribution. The 



bulk of the paper is concentrated in Sec. |V[ where we further analyze the properties of 
the solutions to Eq. (fj). By considering a simpler model of space-independent but variable- 
in-time effective mass we clearly define the exponentially growing "condensate" momentum 
modes. Depending on the number and strength of the condensate modes we obtain different 
distributions for a particle number and distributions over pion species. We emphasize the 
limits when Eq. ([I]) is recovered or when V(f) becomes Gaussian. In the following part of 
Sec. [TV] we consider a full field equation with a time- and space-dependent mass parameter 
and address the exactly solvable case of a mass parameter abruptly changing and returning 
to normal. We show that in this picture the condensate modes are still identifiable and they 
still have a characteristic momentum distribution. Final summary and conclusions are given 
in Sec. [V[ 



II. THE LINEAR SIGMA MODEL 

Below we give a short review of the linear sigma model, point out the origin of Eq. (0) 
in the context of pion fields and deliberate the possible form of the effective mass term. We 
assume the usual linear sigma model Lagrangian of the pion and sigma fields [[18] 

£l.s. = \ M ■ (&>*) + \ (<V) (*V) " \ ((^ 2 + ^) " <) + ea , (3) 

where explicit symmetry breaking is introduced by the parameter e. The sigma field has a 
non-zero vacuum expectation value f n that is set by the Goldberger-Treiman relation and 
is related to the above mentioned symmetry breaking parameter as 

.2 t 2 e 



< = ft ~ 



A/. 
The Lagrangian of Eq. (^) produces the following masses of the pion and sigma mesons 

™ 2 = f, ml = 2\fi + ±. 

Jn Jtt 

Equations of motion for each isospin component of the pion field tt t have the form 

d 2 ir T 
dt 2 



- VV + A (vr 2 + a 2 ) - v 2 ) 7T r = . (4) 



In the mean field approximation non-pionic degrees of freedom, like the sigma field, are 
approximated by their expectation value. Then the dynamics of the pionic field can be 
viewed as a field equation of type (0), with a time- and space-dependent mass parameter 



m 2 cS (x,t) = \((n 2 + a 2 )-vl) . 



As in the application to the chiral condensates, we would like to solve Eq. (^) in the 
formalism of quantum field theory with reasonable expectation values (tt 2 + a 2 ). We do not 
solve the problem self-consistently because m e ff(x, t) , being generated by other degrees of 



freedom, is placed in by hand as a mean field. Before and long after the reaction the fields 
are in their ground state so that (tt t ) = and (a) = f n , which is equivalent to 

m eff( X > t = -°°) = m cff( X > * = +°°) = m l ■ 

During the reaction, the behavior of the effective mass is unknown, and in our case would 
be an input to the model. When chiral symmetry is restored the expectation of the a field 
tends to zero, while the effective squared pion mass is positive. A sudden return of the 
effective potential to its vacuum form (quenching) can strand the a mean field near zero 
with a large negative value of the effective pion mass squared, which in the extreme limit 
would reach m* s = —Xvq ~ —m 2 a /2 . A negative value of m1 a would lead to an exponential 
growth of the low momentum pion modes and long range correlations, i.e. the creation of a 



chiral condensate 19 



III. FROM CLASSICAL TO QUANTUM SOLUTIONS 

The generic form of the Lagrangian density of the field if) that has field terms up to the 
second order in its potential energy part is 

£ = \ (d^) (0"VO - J(x, M- ^m 2 (x, t) ^ . (5) 

To clarify the approach, we consider here a bosonic isoscalar field ip . This does not limit the 
consideration and will be generalized later. It is known that the linear terms J(x, £), often 



called a current or force, can always be removed from consideration p0| . The time- dependent 
current term is responsible for the creation of coherent states with a Poissonian distribution 
of particles. Intensive studies of a DCC produced by linear-type coupling of the pions to the 
disoriented mean field have been performed |I5|| . Unfortunately, many problems cannot be 
easily reduced to this linear approximation. As mentioned in the Introduction, for the chiral 
condensate this would require the Introduction of a symmetry breaking isovector current 
J(x, t), whereas introducing a scalar effective mass (a quadratic term in the Lagrangian) 
preserves intrinsic symmetries. 

There is one important feature of the problem with a quadratic perturbation in the 
Lagrangian. The equations of motion for the field, like Eq. (El), are linear for the Lagrangian 
(0); therefore they are identical for the classical fields and the corresponding operators in 
the Heisenberg picture. Thus, an exact classical solution is related to the solution to the 
quantized version of the problem. The step from classical to quantum treatment will be 
considered next. 



A. Parametric excitation of a harmonic oscillator 

The parametric excitation of a quantal harmonic oscillator has been extensively consid- 
ered in the literature starting from refs. pT|-p5|. This corresponds to our problem in the 



case of no spatial dependence of the effective mass. The presence of translational symmetry 
in this case leads to the conservation of linear momentum. The quantum number of momen- 
tum, k, can be used to label the normal modes. Each mode is just a simple oscillator with 



the time-dependent frequency. The simple Fourier transformation from x to k transforms 
the Hamiltonian, that corresponds to the Lagrangian in Eq. (|5|), into a sum over modes (for 
simplicity we keep one-dimensional notations) 



H=-Jdx (j)(x, tf + V 2 ^(x, t) + m e 2 ff (t)^ 2 (x, *)) = g E (V + k 2 ^\ + m 2 ff (t)^ 2 ) . 

(6) 



With this reduction we arrive at the problem of the independent development of a large 
number of quantum oscillators. With the notation uj 2 (t) = m 2 g(t) + k 2 , we have a classical 
equation of motion for each mode k 

ip k + u&t)il> k = 0. (7) 

After the quantization, the ipk arid the corresponding momenta become operators. Assuming 
a particular normal mode k below we omit this subscript. 

Even for the problem of a simple oscillator, the relation between the classical and quan- 
tum solution is quite subtle; not to mention that classically Eq. (|7|) is analogous to a 
standard problem in quantum mechanics of scattering from a potential — uj 2 that is known 
to have only a limited number of cases with exact analytical solutions. 

Classically, the S'-matrix would be defined if we assume that a solution of Eq. (0) with 
the asymptotics in the remote past corresponding to the frequency u>- , 

^(t)=e" iw -* at i^-oo, (8) 

has evolved with time to a final general form of the solution with the frequency uo + , 



il>(t)= J—(ue- iuJ+t + v*e +iuJ+t ) at t -> +oo . (9) 



It is clear that the complex conjugate to Eq. (||) will develop into a corresponding complex 
conjugate version of Eq. (§). Furthermore, the unitarity, i.e. conservation of probability 
applied to Eq. (0), results in a restriction posed on u and v 

M 2 -M 2 = i- (io) 

In order to analyze the quantum version of the problem we use the language of secondary 
quantization, introducing creation and annihilation operators for a single harmonic oscillator. 
The time dependence of a quantized field coordinate ip(t) in the Heisenberg representation 

as t — > — oo is 

m = -J=( oe ---< + aV^), (11) 

where the operators a and a* do not depend on time and define the "in" state of the system. 
We will define the "out" state using operators b and tf as 



V>(t) = -==(fee-^ +t + bU iu+t ) at t -> +oo. (12) 



Being supported by the fact that the field ip(t), even in its quantum form, satisfies Eq. ((7|), 
and with our assumptions about the classical solution (^ [|) we can relate operators a and 
a) with b and tf , 

b = (ua + va^ , b ] = (u* a f + v* a) , (13) 



with the same parameters u and v as in the classical solution. As seen ( p"3"D the "in" and "out" 
states are related by a Bogoliubov transformation. The condition that the commutation 
relation 

[6,fet] = (|«| 2 - M 2 ) [a, 0+] = 1 

is preserved coincides with Eq. (|ToD . 

It is possible to build an S- matrix corresponding to the transformation in Eq. (|13|). 
however the solution is very technical, and of limited use. We can use Eq. ([13|) directly to 
answer simple but relevant questions. The probability amplitude to have n final quanta if 
there was vacuum in the initial state is given by 

C n = a (n\S\0) a = a {n\0) b , 

where the S-matrix is a transition matrix from the "in" state to the "out" state. Then 

oo 

|0>6 = E C n\n)a, 
n=0 

and recursion relation for the coefficients C„ 



C n+ 2---\ r ^C n (14) 

u V n + 2 

may be found from the definition of the vacuum and Eq. (|i~3|): 

oo 

&|0>6 = Y,°n (u Vn \n - l) a + v Vn + 1 \n + l) a ) = 0. 

ra=0 

Finally, normalized to unity, the transition probability from a vacuum to a state with 2n 
bosons is expressed by 

l^| 2 = ^) = ^|^p"v/W, (15) 

where the parameter p = \v/u\ 2 represents the reflection probability for Eq. (0) if it is viewed 
as a Schrodinger equation for scattering at zero energy off a potential — u 2 . Conservation 
of probability, Eq. (0), limits the values of p to the region p G [0, 1) . The distribution of 
particles produced from vacuum Eq. (|15|) peaks at zero and the average number of particles 
produced is 

00 o 

n = E 2nP(2n) = — — . (16) 



n=0 



1-p 



It is seen from the above expression that the number of particles diverges as p approaches 
one. This will be a region of interest in our later studies of chiral condensates. 
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A rigorous and detailed discussion of all the features in the transition from classical 
solutions of Eq. ([7]) to a problem of a quantum oscillator with a time-dependent frequency 
can be found in pl|-p3|. The corresponding Schrodinger equation for the wave function of 



the quantum oscillator is (h = 1) 

,d$(x,t) 1 d 2 <5>(x,t) 1 2 . . „ , . 

— 5T i = -2-^ Li + 2^ ( «^* (l -' ) ' (17) 

where coordinate x originates from the field "coordinate" ip, Eq. (0). Let some function 
ip(t) be a solution to classical Eq. (0) with c<j(t = — oo) = cj_ and cj(t = +oo) = u> + . 
This function ip(t) can be written in terms of the real part r(t) and phase j(t) so that 
i/)(t) = r(t) e* 7 ^ . With a direct substitution it can be shown that 



$0M) = -r=exp [i— )x(y,r) (18) 




is a solution to the quantum equation (|rT|). Here a rescaled length y and time r are in- 
troduced as y = x/r(t) , r = j(t)/u>-. The function xO^jO * s a standard solution of the 
Schrodinger equation for the oscillator with a constant frequency uj- , 

• dx(y, t) 1 d 2 x(y,r) 1 2 

^^^ = -2^^ + 2^ ^' T) - (19) 

The initial conditions for the solution x{Vi T ) should be set at t — ► — oo by the known 
conditions on $(x, t) . Equation (|15D can be generalized utilizing this more general technique. 
It can be shown that the probability for a transition from the state with m quanta into the 
final state of n quanta is given by 



mm(m, n) 
maxfm, n) 



1-p 



P (Z+n)/2^\J l ~ P) 



(20) 



where Pjf are the associated Legendre polynomials and n and m have to be of the same 
parity. The established bridge between classical and quantum solutions is a remarkable 
achievement, but unfortunately analysis of Eq. (|I8"D does not present a pleasant task. 



B. Infinite number of mixed modes 



1. General canonical transformation 



The situation becomes more complicated if we return to the general field equation given 
by the Lagrangian density of Eq. (|), and assume that the modes cannot be separated. We 
keep assuming J(x, t) = which is relevant to our particular problem but in general this 
factor can be included back in the discussion without much complication. The problem of the 
system of coupled oscillators has been discussed in great detail in [|4]]. The transformation 
analogous to Eq. ([13|) now takes an iV- dimensional symplectic form, where N is the number 
of coupled modes 



N N 



h = J2 ( u kk>ak> +ffefe'4') > & l = 5Z ( M fcfc' a l + v lk' a k') , (21) 



fc'=i fc'=i 



It is convenient to regard the operators a k as components of an N- dimensional vector, and 
the numbers Ukk 1 and Vkk> a s N x iV matrices. Eq. (511) in matrix notation takes the form 



b =S~ 1 aS = ua + vat, V = S' 1 a^ S = u* a 1 + v* a . (22) 

Using the properties of the matrices u and v shown below, it is straightforward to check 
that the inverse transformation is given by 

a =SbS- 1 = uH - v T b\ a) = Stf S~ l = u T b ] -vH, (23) 

where u*, u^ and u T have their usual meanings of complex conjugate, hermitian conjugate 
and transpose matrices, respectively. Further, we will also use an inversion denoted as u^ 1 . 
Similarly to the one- mode example of Eq. (|10D , matrices u and v are subject to conditions 
that arise from the fact that the commutation relations have to be preserved. It follows from 
Eq. <^j that 



u u' — v v ' = 1 , u v T — v u T = , (24) 

and from the inverse transformation, Eq. ( |23| ) 

u ] u - v T v* = 1, u t v-v T u* = 0. (25) 

2. Transitions between coherent states 

In the case of many mixed modes it is still possible to obtain recursion relations similar 
to the single-mode situation of Eq. (Ill) but they become difficult to analyze. Next we 



will discuss the approach found in |E3| with a different technique that uses the transitions 
between coherent states. We define a single-mode coherent state \a) in the usual way as 

oo n 

l«>=e-'"' 2 / 2 £^>- (26) 

In mathematics sums as in Eq. (|26| ) are often called generating functions. It is convenient 



to introduce states with a different normalization [26] as 

\\a) = e |a|2/2 |«) = e aat |0). (27) 

When applied to the states fl27|) , the creation operation is equivalent to the derivative, 

d 
aU\a) = —— I la) . (28) 

da 

We continue to use the notation |ot) for a multidimensional coherent state, which should 
be interpreted as a product of the single-mode states, a. is a vector and the derivative d/dct 
is understood as a gradient in iV-dimensional space of modes. We are interested in finding 

8 



the matrix elements of the evolution matrix S that implements the unitary transformation 
between initial and final states of the system. The evolution of the initial state is quite 
complicated and unless v = the coherent state does not stay coherent as in general it 
evolves into a so-called "squeezed" state [p7| , p8[] . Nevertheless it is possible to obtain an 



analytic expression for the matrix elements of the evolution operator S between coherent 
states. To proceed in this direction we will consider the action of the matrix S on the 
creation and annihilation operators given by Eqs. ( |2"2"| , |2lf ), that actually serve here as the 
definition of the evolution matrix. By acting on the complex conjugate form of the first 
equation in Eq. (E^) with (f3 \\S from the left and with ||tt) from the right, and utilizing 
Eq. fl28[ ) we arrive at the following differential equation 

(p* ~ u*^- - v*a) ((3\\S\\a) = 0. (29) 

In a similar manner from the first equation of Eq. (|23| ) we obtain 

U - u^ + v T pA (P\\S\\a) = 0. (30) 

The solution to the differential equations (|29|,|30D determines the transition amplitude up to 



a normalization constant C(u,v): 

(P\\S\\a) = C(u, v) exp (tf (u^ 1 a + ~ /3 f v (u*)' 1 P* - )-cx T (u*)' 1 v* a\ . (31) 

This solution can be checked directly by substitution into Eqs. ( |29|J3"0| ) and utilizing the 
observation which follows from Eqs. ( P^ , P3|) that matrices v (u*) _1 and (u*)~ l v* are both 
symmetric. The normalization constant may be obtained using the completeness of coherent 

states 



-1/4 



C{u,v) = (det(W)) . (32) 

The evolution matrix given in the form of Eq. ( |3"T| ) can be transformed via a Taylor expansion 
in the particle number basis using the definitions in Eqs. 



_ N n*n k N n t 

(P\\S\\a) = £ (n*|5|n{> U^U^T. ( 33 ) 



{n k }{n\} k Vn k - i yjn\\ 

In general it is quite complicated to give a finite expression for the coefficients of the Taylor 
series arising from the multi-variable Gaussian. Expansions of two-variable Gaussians are 
known to be of the form of Legendre polynomials which give rise to Eq. (p0|) . Nevertheless, 
the algorithm for the expansion is straightforward, first the exponent should be expanded 
in terms of its argument and then each term can be expanded into a final sum with a 
generalized binomial expansion. 

By considering a transition from the vacuum state we set all terms with ex to zero in 
Eq. (|31"|). In this case everything is completely determined by the matrix v (w*) -1 . For the 
question of multiplicity distributions in the one-mode case, we note that the relative phase 
between u and v was of no importance, a single parameter p determined everything. We 
will further see that a similar picture holds in the general case, and only one matrix vv* 
is needed to find the particle distributions. With the help of Eqs. (JM, |23|) it can be seen 
that the matrices uu^ and w* can be diagonalized simultaneously and the eigenvalues of 
vv ^(W)'" 1 ' = v (n*)~ 1 (t> (u*)^ 1 )^ form a set of parameters p different for each mode. 



3. Multiplicity distributions 



Despite the complicated form of a general expression (^T|) one may calculate the mo- 
ments of multiplicity distributions in a straightforward manner. Let the initial state \i) be 
characterized by a diagonal density matrix n°, a set of numbers of quanta {n%} in each mode 
k on the diagonal. The average number of quanta in a final state |/) is determined by 

nk = (f\ala k \f) = {i\b\b k \i) . 
With the above assumptions and Eqs. (E3, 1231), the density matrix n of a final state is 



h = un°u* + vn°v* + vv* . (34) 

Throughout the rest of the work we will concentrate on the situation of a particular interest 



when particles are created from the vacuum and the first two terms in Eq. (34) are identically 
zero. The average total number of particles in this case may be expressed in a simple matrix 
form 

rkotax = Yl n k = ^(v*v) , (35) 



which is consistent with Eq. ([Tq). Higher moments of the particle distributions can be 
calculated in the same manner. Unfortunately, the calculation of an arbitrary moment 
requires path integration techniques while using Wick's theorem in the normal ordering of 
operators. Low order moments can be directly calculated, for example 



n 2 t otai = n t 2 ota i + 2n to tai + 2 Tr [(vv* ) 2 ) , (36) 

^total = <taJ + 6<tal + ^total + 6n total Tr ((w* )2) + UTl (( w t )2) + gTr 



'v v* ) 3 



(37) 



Equation (|36|) can be identified as a super- Poissonian distribution of particle pairs. It follows 
from the above expressions that particle production is determined by the matrix vv* which 
is related in a simple way to uu* , Eq. fllOf) . As was mentioned before, the hermitian matrix 
vv* can always be diagonalized with diagonal elements being average numbers of particles 
n^ in each eigenmode q. This diagonalization allows us to view particle production as a 
production from independent modes. A connection can be established to the one-oscillator 
case, previously considered, by defining a set of parameters p as p q = 7^/(1 + n^) , that are 
the eigenvalues of vv*uu* . Considering the number of particles produced in each eigenmode 
we can restrict ourselves by the modes with p q — > 1 , that dominate the particle production. 
We will later refer to these modes as condensate modes. In the following section, it will be 
shown that the distinct physical feature of these modes is that they produce an exponentially 
large number of particles. 

Next we consider a number of condensate modes with equal parameters p q . For one 
mode the answer is in Eq. flTo]). For the case with several condensate modes one needs to 
know the distribution for the total sum of the particles, which is given by the convolution of 
the corresponding probabilities. As for the Fourier transformation, a convolution of several 

10 



distributions results in a product of their generating functions. The generating function for 



the distribution (fj) is 

n=o v L -py 

We obtain total probabilities of particle production for any combination of species or any 
number of modes convoluted together, as a Fourier expansion of a product of corresponding 
generating functions. Suppose we are looking for a distribution P/(2n)that gives a probability 
to observe 2n particles appearing in I single-mode distributions ([HJ). The Taylor expansion 
of (1 — py)~ 1 / 2 gives 

P( , n) (* + 2(n-l))!! l/2 

Pl{2n) ~ 2"n!(l-2)!! P (1 ~ p) ' (38) 

where the double factorial should be understood appropriately for odd and even I . 

According to the central limit theorem, only the distributions corresponding to a small 
number of condensate modes have a shape that is very different from a Gaussian, see Figure 
HI below. If parameters {p q } for the condensate modes are very different then only the 
important modes that produce many particles can be isolated reducing the problem back 
to the case of several almost degenerate modes. The method of generating functions is 
also convenient in discussions of distributions over species. For example, the distribution 
for the total number of pions, regardless of the isospin projection, can be obtained as a 
distribution of the sum of three species, via convolution. Alternatively, one can think of 
the total condensate modes for pions as just a sum of numbers of condensate modes for all 
species. 

As an example we apply Eq. (|38| ) for one mode and three pion species. The probability 
of having 2n neutral pions is 



(2n)\ 
2 2n (nV 



P%2n)=P 1 (2n) = -^-(l-p)^p\ (39) 



the probability for observing 2n charged pions is 

P+(n) = P-{n) = P chargcd (2n) = P 2 {2n) = (1 - p)p n , (40) 

and finally the probability of having total 2n pions is 

P-(2n) = P 3 (2n) = ^±1^ (1 _ p? l* f . (41) 

We note the convolution of three distributions in Eq. ( p9|) that give rise to Eq. ( pEip produces 
a peaked curve with maximum at particle number 

3p-2 

^max ~Z ■ 

1-p 



The probability distributions ([39], [ID], ^T[) are shown in Fig. |T] with solid line, long dashed line 
and short dashed line, respectively. The figure displays a critical situation when parameter p 
is 0.999. This corresponds closely to the first mode in the example of a square perturbation 
considered in the next section. 
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FIG. 1. The particle number distribution for neutral pions (solid line), charged pions (long 
dashes), and for all pions (short dashed line). All curves are normalized to unity; note that the 
total number of pions and number of neutral pions both are always even, while any integer number 
of charged pions appearing in pairs is allowed. The parameter p = 0.999 . 

Now we can return to original questions: What is the probability to observe a certain 
fraction of neutral pions given a total number of pions detected? Is it possible to use this 
distribution for the detection of chiral condensate? In the one-mode approximation, given 
a total number of pions produced 2n t , the normalized probability of observing 2n neutral 
particles is 



V 



/=- 



P°(2n) P ch ^ ed (2n t - 2n) (2n) 



P tot (2n t 



(2n t )\ 




^1 2 



2nt—n 



(42) 



Finally, if both n and n t are large, which is a good approximation in almost all regions, 
Stirling's formula may be used, and one finds 



W) 



(2n t + 1)77 



(43) 



which coincides with Eq. ([I]) if the normalization over n t + 1 discrete points between zero 
and one is switched to the integral over / e [0, 1]. 

Unfortunately, this result does not hold for all situations when several modes are par- 
ticipating together. The single-mode result fl4~3"D may be invalid when the largest eigenvalue 
of the matrix vv< is exactly or almost degenerate. Physically, this may be due to some 
symmetry for example. One- and two-mode cases have their probability peaked at zero, this 
is no longer true for the convolution of three or more modes, Fig. [I]. Eq. (^) can be applied 
to give a result for any number of modes j participating in the condensate, assuming they 
all have equal strength p . The distribution of neutral pions is given by 
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T>j\f 



n 

m 



P j (2n)P 2j (2n t -2n) 



With a Stirling's formula and the assumption n ^> j , Eq. 

(3J/2-1)! 



(44) 



vAf) 



(j/2-1)!(j-1)! 



|) can be simplified to 
p/2-i (1 _ f y-i . (45) 



For comparison, the probability distributions for one, two, and three modes in the condensate 
are shown in Fig. FJ. 




0.2 0.4 0.6 0.8 

neutral pion ratio f 



FIG. 2. The probability V(f) that a given neutral pion fraction / is observed. The three curves 
display cases of one condensate mode (solid line) , two energy-degenerate modes (dotted line) , and 
three modes (dashed line). 



The distributions of Eq. (|45f ) have maxima, average and widths as follows 



/«. 



3J-4 



3 /=l/3, ; 2 -/ 2 



27 j + 18 ' 



where the first equation can be applied with the restriction j > 2 . It follows from the central 
limit theorem that for a large number of modes one should expect a Gaussian distribution 
that tends to a ^-function. Thus, as the number of modes in the condensate grows, there is 
a fast transition from a 1/yff type behavior to a sharp peak at / = 1/3 . 



4- Dynamics of matrices u and v 

As some final remarks about the link between classical and quantum solutions to the 
wave equation arising from the Lagrangian in Eq. (13), we would like to discuss the technical 
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question of constructing matrices u and v . Unfortunately, Eq. (|2|) has no general analytic 
solution even classically. The Green's function formalism reduces the problem to a Fredholm 
type integral equation. The exact solution can be obtained only in special cases, for a 
separable kernel p0| , or in the sudden perturbation limit to be considered in the next 



section. However, numerical studies of Eq. (fj) seem to offer a great chance of success. 

It would be good to obtain equations for u and v that are still exact but written in the 
form convenient for numerical work. Let us quantize the field at every intermediate stage 
with bare particles PU| so that 



i}{x,t) =£ —1= (&,(£) e^"^ 4 + bl(t)e- ik * +i " kt ) , (46) 

k \J luj k L K y 

where b k (t) is now a time-dependent annihilation operator. Eq. ([SJ) is written in an 
interaction representation, the explicit time dependence of the free field is given by the 
exponents, whereas creation and annihilation operators absorb the remaining nontrivial 

1 /9 

time dependence. In Eq. (|46| ) the variable uj k = (k 2 + rn 2 e ff(t — > — oo)J is the initial 
time-independent frequency, thus there is no problem if at any point in time the effective 
mass goes through zero. As time goes to infinity i/j(x,t) becomes the final "out" state, with 
operators b and tf defined as before. For one-dimensional case L is the quantization length. 
For further simplicity we to denote 

m 2 eff {x,t) = m 2 eff {t^ -oo) + n(x,t), (47) 

this allows the separation of the interaction Hamiltonian. The perturbed Hamiltonian from 
Eq. (H) expressed in terms of b(t) and b^(t) is 



H**(t) =£ (n kV (t) (bi(t)b k ,(t) + 6 -f-) + lA k At)bl(t)bUt) + ~A* kk ,(t)b k (t)b k ,(t) 



(48) 
where the matrices Q and A are determined as follows: 

n kk ,(t) = or \ a (u(k - k',^'-^*) , A kk ,(t) = n (- fc, -^_*) e *v+«*)t. (49) 

2LyJu k L) kl v / 2Ly/uJ k uj k/ 

In the above expression U(k, t) is a Fourier image of II(a;, t) , determined as 

U(k,t) = Ii*(-k,t) = J U(x,t)e ikx dx, 

and 3? denotes a real part of the expression. Utilizing the Hamiltonian equation of motion 
in the interaction picture 

i^b(t) = [b(t),H mt -(t)} = nb(t) + Atf(t) 

at 
it is possible to show that if b(t) is defined through the initial operators a and a^ as 

b(t) = u(t)a + v{t)a\ (50) 
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where u and v must satisfy the matrix equations 

i—u(t) = n(t)u(t) + A{t)v*{t), (51) 

CiiL 

i j- v(t) = Q(t) v(t) + A{t)u*{t) . (52) 

(I I 

At infinitely large times when the perturbation II goes to zero, the right hand sides of Eqs. 
( |5*T| , |5*2] ) vanish and u and v become time-independent. In order to obtain matrices u and v 
for Eq. (P2"D one has to solve the above equations with initial conditions u = 1 and v = . 



IV. APPLICATION TO CHIRAL CONDENSATE 



A. Separable modes, space-independent effective mass 

To illustrate the machinery developed in the previous section we start with simple cases 
when the classical solution is known analytically. The first example is the space-independent 
field II(x, t) = Il(t) , see Eq. ([47D, i.e. the situation when the perturbation is uniform in a 
box to which the entire pion field is confined. The wave vector k is a good quantum number. 
Particles get produced independently in each mode labeled with k, and production is deter- 
mined by the classical reflection probability p^ = |fk| 2 /|wk| 2 j Eq. fll5|). The distribution of 



particles for a single mode, Eq. (15), is a decaying function that has a maximum at zero. 



Its behavior can be approximated with Stirling's formula as 



P(2n) ~ ^ > 



where P(2n) denotes the probability of creating 2n particles from the vacuum. The distri- 
bution in Eq. (|15|) for a value of p = 0.999 is shown in Fig. [I] as a solid line; some additional 
discussion is given below. The average number of particles produced in Eq. (|1^) is generally 
quite small unless we are in the condensate region when p — > 1 . With the assumption of 
independent modes Eq. (|) reduces to 

- d 2 7i T /dt 2 - U(t)n T = (ml + k 2 )n T . (53) 

Eq. ( [5*3"D is written in a form of the Schrodinger equation for scattering from a potential 

barrier of height —II at "energy" m^ + k 2 . The parameter p we are looking for, which links 

classical and quantum pictures, is the reflection coefficient for this scattering process. 

Let the effective pion mass change and then return back to normal in a step function 

manner. This corresponds to the scattering off a rectangular potential barrier. The situation 

when the tunneling is involved, is relevant to the case of a low momentum mode being 

amplified as the index of reflection is rapidly increasing. We assume that the perturbation 

Il(t) has a non-zero value II only for the time interval t G [0, T] . The reflection probability 

for this scattering potential is 

Kl 2 
P^ ~ 7T1 — 12 ' 

1 + \Vk\ 2 
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and the average number n-^ of particles created in the mode k is given by 



n 2 



n k = Vk 



A(ml + k 2 )|m 2 + k 2 + III 



sin TJml + k 2 + II 



(54) 



Another form of the perturbation II (t) for which an exact analytical solution exists is 
the Eckart potential p9j 

n 

{t) ~ cosh 2 (t/T) ' 

where II = 11(0) is the minimum value of U.(t) and T is the time scale of perturbation. The 
resulting form of the reflection probability is 



Pk 



1 + cos (vr V4ITT 2 + 1 



cosh f 2 7r T 



mi 



cos (it V4:UT 2 + 1 



Both forms of the perturbation H(t) with parameters that we use below for our numerical 
estimates are plotted in Fig. |3], right side shows the Eckart potential, rectangular barrier is 
on the left. 
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FIG. 3. Eckart (|II| = 1 (GeV) 2 , T = 0.5 fm/c) and rectangular (|H| = 1 (GeV) 2 , T = 1 fm/c) 
perturbations are shown on the right and left panels, respectively. 

Fig. f| shows the average number of particles produced in the two models described 
above. For our estimate we made the following choice of parameters. The effective mass 
drops to the value of m 2 eS = — m 2 /2, correspondingly the parameter II is chosen at —1 GeV 
for both models. The time scale, given by parameter T, is 1 fm/c and 0.5 fm/c for the square 
barrier and the Eckart potential, respectively. It is important to notice that even though 
the graph is plotted over a continuous variable k , our finite size spatial box allows only 
discrete values of momentum in each direction. With the expected interaction region size 
of 2 fm, the lowest momentum in a cubic well of size L is |k| w \^3ii/L, numerically this 
is about 540 MeV. According to Fig. [|, both models predict around several hundreds of 
particles. The next higher-lying modes, have significantly smaller numbers of particles. The 
picture presented shows that one can only hope to have very few modes that actually form 
the condensate, as the number of particles falls drastically for higher momenta. This is 
consistent with the argument that in order to get a noticeable condensate one should have 
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energy of the mode dipping below zero. The square barrier model provides a simple estimate 
for the number of particles produced if the mode just touches zero, 



ITT 2 



n 



mlT 2 



6. 



where the sigma mass m a is taken 1.4 GeV and the interaction time T is around 1 fm/c . 
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FIG. 4. The average number of particles produced in the mode k as a function of k for square 
potential (solid line) and Eckart potential (dashed line). The parameters of the perturbation are 
chosen so that II = 1 GeV 2 for both models, and T is 0.5 fm/c and 1.0 fm/c for Eckart and square 
potentials, respectively. 

For more complicated perturbations it is possible to use the WKB approximation in 
order to determine the reflection coefficient of classical Eq. (|7|) for tunneling. With the 
notation 

f = 2 exp ( j°° \u(t)\ (-uJ 2 (t)) dt\ , 

where the Heaviside theta- function limits the integration to the region of negative w 2 (t) , we 
have the average number of particles given in the semiclassical limit of £ ^> 1 by 



_ £ 2 



(55) 



Concerning the total particle distribution one has to add all particles from all modes and 
consider the distribution as a superposition. Based on the results of the above examples 
and Sec. 2 we can conclude the following. About a thousand neutral pions with momentum 
540 MeV/c are produced by the lowest mode and are distributed according to Pi(2n), Eq. 
(J38[). The distribution over species is V\ ~ 1/y/J , Eq. (|45|). This is the single-mode 
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result described by Eqs. (|39| , |40| , f41|) and shown in Fig. |1|. Geometry is crucial here, as 
it determines the energies and degeneracies of other modes that may or may not compete 
with the lowest mode(s). Higher momentum modes in the box do not "condense" if, in 
the scattering picture, their energy is higher than the barrier and reflection is negligible. 
These modes, even jointly, may on average produce just several particles. The case of 
the non-condensate particle production will be considered with a better model in the next 
subsection. 

These simple examples are still far from realistic. One of the major failures is that pions 
from the square box are not real pions and therefore all excitations that we obtain need 
to be projected onto final pion states given by the plane waves of the entire space. This 
projection will produce momentum spread for particles, that will carry the characteristics 
of each condensate mode. Nevertheless, these models produced reasonable results and what 
is more important they have identified the physics of the process. 



B. Time- and space-dependent perturbation 



1. General solution for a step-like temporal perturbation 

Here we will solve the perturbed Klein-Gordon equation for the pions in a more realistic 
case, where the effective pion mass is both space- and time-dependent, 

_1 _ v 2 tt + {ml + n(x, tj) n = . (56) 

Despite the fact that this dependence is put into the model by hand, the resulting features 
can be quite general. In order to keep our solutions analytic we choose the perturbation 
n(x, t) as a step function in time 



n(x,t) 



n(x) for te [0,7] 
otherwise 



This choice will allow us to solve Eq. (EH) classically and to construct matrices u and v that 
control the evolution of a plane wave. For most of the discussion the form of n(x) is left as 
general, but in the last part for the numerical results we take it as a spherical square well. 
The wave function tt(x, t) can be found in each of the time regions, and the solutions 
should be smoothly matched keeping 7r(x, t) and the derivative dir(x.,t)/dt continuous. In- 
troducing a separation of variables as 7r(x, t) = X(x) T(t) we obtain equations of motion 
for X and T in all three time regions. In the perturbed region t G [0,T] , Eq. ( J56f) becomes 

- v 2 ;t(x) + n(x) x(x) = s x(x) , (57) 

d 2 T(t) 



dt 2 



W 2 T(t) = 0, 



where the dispersion relation is 



W 2 = m 2 + £ . (58) 



The unperturbed form of the free space wave equations at t < or t > T is 



v 2 x°u) = \k\ 2 x°M 



(59) 



d 2 T°(t) 
dt 2 



with the usual relation 



UJZ 



u£T{t) =0 



m^ + Ikl 



Eq. ([59"D has simple plane wave eigensolutions that we denote as |k) with positive eigenvalues 
|k| 2 . The perturbed Eq. (57) with negative II may have negative energy bound states as 
well as the usual continuum states with positive energies, see Fig. ||. For such a bound 
state, W 2 , Eq. (|5|), becomes negative and W = iQ is imaginary. 



n(x)A 




FIG. 5. The schematic representation of the perturbation II(x), with one condensate bound 
state. 

We denote the eigenstates of Eq. ( |5TD as \n) and corresponding energy as £(k) . Here n 
is a set of quantum numbers labeling the eigenstates. We also assume that both |k) and \k) 
are properly normalized and form complete sets, 



<k|k'> = <5 k)k , 



[K\K 



s K , K/ , i = El k >< k l = EI«)(*l- 

k k 



To determine the Bogoliubov transformation matrices u and v we consider the evolution of 
the wave function 



|7r(t)> = e 
through an intermediate stage t e [0, T] 



-iuj^t 



k) , for t < 



t(*)> = E K.e-^ + 6 kK 



iW K t 
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into a final state 



\<t)) = E \ — («kkfe-*^» + vwe^'t) Ik'). 



Continuity of the wave function and its derivative with respect to time allows the determi- 
nation of unknown coefficients ak K and b^ « as well as the matrix elements of interest Vkk' 
and «kk' • The result is, both for real and imaginary W K , 



iui^T 



^kk' 



u; k ' e * 



U} k 



E< k V)W 



l + ^)cos(^T)-,(^ + ^)su l( nV/M 



(60) 



and 



^kk' 



' > (kkj(ftk) x 



i-^)oos(w B r)-i(^-^)si 11( ir,:n 

^k'/ \W K Wfc ' 



(61) 



In order to see whether the condensate was formed we first address the question of the 
number of produced pions. Possible bound states in the solutions of Eq. fl57|) should 



be carefully treated. As seen below, bound states in Eq. fl57|) with the energy £ below 
— m\ produce an exponential growth of the particle number with time which is our starting 
criterion for the search of the condensate. The average number of particles may be then 
expressed using Eq. (|35f ) as 



n 



E l^k 



k,k' 



•"I 2 = ^ E i A «« cos ( W - T ) cos ( W *' T ) + B **t sin ( W - T ) sin (W$T)} , (62) 

re , re' 



A KK / - 2 [I KK , I KK , - 8 KK * 



r(i) 2 



"k w.' 



w.wj W w« 



w " ^'U.+4-vV^i, 



where we assumed that the states are chosen in such a way that the amplitudes (k|/c) are 
real, and we introduced the notations 



r(i) 



^re'=^re = EWkK(k|/,') 



and 



j(-i) = j(-i) 



$>|k>-<k|«/). 

k ^k 



Depending on a particular level, P^ may be real or imaginary if in Eq. (|53) the energy £ is 
greater or less than m 2 ; nevertheless Eq. ( |62"D works for both cases. 

In practice it is convenient to express the overlap of (k|«) using Eqs. (0, ^) for the 
states |«) and |k) . Multiplying Eq. ( p7|) by A" and Eq. fl5D| ) by X and subtracting results 
we obtain 
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(£ - k 2 )X X° = U(x.)X X° + V(XVX° - X°VX) . (63) 

Integration of Eq. (^3|) over all space produces a useful result 



(k|«) = — —J d i xU(x)XX . (64) 

For the continuous spectrum, this expression contains an additional term with 5(8 — k 2 ) . 

Due to large oscillations of the trigonometric factors in Eqs. (|6~lj ) and (|62l) , one should 
expect considerable particle production only in situations with imaginary W . This observa- 
tion makes it natural to separate the sum in Eq. fl6~2|) into several contributions depending 
on the intermediate state \k) . There are exponentially rising terms that involve transitions 
through the states with imaginary W; it is important that the bound states are discrete 
and their number is finite. Another contribution is of all intermediate states that lie in the 
continuum. This second contribution is always present for any perturbation even with no 
bound states. With a non-zero pion mass there maybe a number of discrete states that 
have real W but these states have negligible contributions compared to the states from the 
continuum. In the remaining part of this subsection we will mainly draw attention to the 
first two cases which we call "condensate" and "non-condensate" pion production. 



2. Condensate pion production. 

In the following picture we assume that there is one state \kq) with negative energy so 
that W K0 = ifl , and in the summation in Eq. (|61~1) the term that involves this condensate 
state \kq) is dominant. Thus, the distribution of particles and all other properties are the 
same as in the single-mode example, i.e. as for the parametrically excited single oscillator. 
According to Eqs. Q35|) , (p6|), (J37p and (p4]) , the traces of w^ and other higher powers 
of this matrix completely determine all moments of the particle distribution and thus the 
distribution itself. Therefore, it is sufficient to consider the eigenvalues of a hermitian matrix 

Under the assumption that only one term in the sum in Eq. ([H]) is important, we can 
express the matrix element of vv* as 



{W}k,k' = a - +by/UkU}K + cJ— + c* J—, (65) 

V^k WU V UV V ^k 

where the introduced coefficients are 

a = - (/« cosh 2 (ftT) + tfl^ sinh 2 (ftT)) , (66) 

b = - (/ ( - 1} cosh 2 (fiT) + /« sinh 2 (fiT)/fi 2 ) , (67) 

c =i(i[j( 1 Vn + / ( - 1) fi] -l) , (68) 



and 



I'-^EKkiKo)! 2 -, i (1) = E l(k|«o)l 2 ^ 
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We wrote an element of the matrix vv^ in the form of Eq. ( J65|) in order to show that it 
consists of only four factorizable terms. In accordance with a general theory of factorizable 
kernels, the only non-zero eigenvalues of matrix vv* are the eigenvalues of the following 4x4 
matrix 



al { -^ b c 6/C- 1 ) 

a &/« clV) c * 

al^ b c bl(~V 

a 6/W c/W c* 



\ 



(69) 



Having the second two rows the same as the first two in the matrix (|69"D reduces the actual 
number of non-zero eigenvalues to two, and the secular equation for the eigenvalues A 



A 2 -ATr( W t) + l(/(i)/(-i)_i)2 
16 







(70) 



The trace of the matrix vv* in this expression is the same as the average number of particles 
produced by one condensate state |«o) , 



n = Tr(tnr) 



cosh 2 (fiT) 



sinh 2 (^T) f J( _i)2 fi2 



1 
2 



(71) 



The result (|7T|) can be also obtained directly from Eq. 
In the limit of exp(2f2T) ~^> 1, Eq. (|7l"D becomes 



n 



n 2QT 
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i^n + 






(72) 



This approximation is relevant to our problem as we wish to determine the exponentially 
growing condensate modes and therefore we choose a physical environment for which this 
is true. The limitation exp(2QT) ^ 1 is then a criteria for the environment of chiral 
condensate. As seen from Eq. (|70|) , in the same limit one eigenvalue of matrix vv' becomes 
exponentially large and the other one goes to zero. This is important because with only one 
nonzero eigenvalue we recover the single-mode situation, discussed in the previous sections 
which leads to the distribution of particles of the form of Eq. (J39|) and a 1/y/J form for the 
distribution over the pion species. 

As a concluding statement we stress here that in principle just making the approximation 
exp(2QT) ^> 1 allows us to take one term in the sum of Eqs. fl60[ |B"l"| ) from the very 
beginning, and the same approximation was used again in the end to purify the single 
mode. Therefore the condition exp(2f2T) ^> 1 must be a clear indication of a condensate in 
the bound mode with relative energy Q . If there are several bound condensate states with 
close energies Q we would again recover the system of several modes that was discussed 
before. 
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3. N on- condensate pion production 

To make the picture complete we have to estimate the number of particles that are 
produced from all the modes not involved in the condensate. The low-lying negative energy 
states that have energies above the pion mass do not make a considerable contribution. They 
may be estimated by the same method that was used for a condensate states just having 
oscillating exponents instead of exponential growth in Eq. (|72|), 



If,, JW\ 

" ~ h r i)n + V) • < 73 > 

This source does not produce many pions, as will be seen from the numerical example in 
the next section. 

The contribution to the number of particles from the continuum presents a greater in- 
terest. First, as the number of states is infinite, unlike the previous case of bound states, we 
may have a significant contribution. Secondly, it is an important practical question because 
the continuum states always produce pions even if the condensate is absent. It follows from 
Sec. |l I| that the charge distribution of pions from many modes of almost equal strength 
is Gaussian. To perform a particle number estimate we again address Eq. ([7ip. It is rea- 
sonable to assume that oscillating terms corresponding to different arguments W K and W K > 
average out to zero and the contributing terms are those that are in phase, corresponding 
to the same k . For these terms we take sm 2 (WT) rs cos 2 (WT) rs 1/2. Applying these 
approximations to Eq. (ffll) we obtain 

HEKkwr^^- (74) 

1 k k w k yy K 



As seen from the structure of this equation, the arbitrary ^-contribution in Eq. (pty does 
not influence the result for n . 



C. Spatial spherical square well 

The spherical square well potential of a finite depth is the simplest spatial perturbation 
n(x) that has an analytical solution. We assume that 

-V for |x| < R 

LI(x) \ 

otherwise 

The rotational symmetry makes angular momentum a good quantum number, and the 
remaining radial dependence can be expressed in terms of spherical Bessel functions. Due 
to the largest exponential enhancement the deepest level is expected to produce most of the 
contribution. Thus, we concentrate our attention on an s-wave bound ground state that 
will be a dominant condensate state in the process. By the virtue of symmetry, all pions 
produced from this state will have a spherically symmetrical s-wave spatial distribution. 
Instead of plane waves it is convenient to quantize spherical waves in a large sphere of 
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radius L . This leads to the substitution of the former basis |k) by \k,l,m) where k = |k| 
and I and m are the orbital momentum and its projection, respectively. Needed overlaps 
between perturbed and non-perturbed states can be computed using the one-dimensional 
version of Eq. (|64] ) where the right hand side is obtained by integrating Eq. (^) from zero 
up to the size R giving 

V (X(r)d/drX°(r)-X°(r)d/drX(r)) 



(k|K) " £-W E-k^ + V 

In particular, for a bound (£ < 0) s-state |k ) Eq. ( |75]) gives 

2a 2 acot(aR) sin(kR) — kcos(kR) 



--H 



(«o|fc,0,0) 



a 



L sin(ai?) (k 2 



a^ 



(k 2 + a 2 cot 2 (aR)) VaR- tan(ai?) 



(75) 



(76) 



where a = JV —\S\. This overlap is normalized to one as a sum over all momenta k = 
7fn/L , where n is a nonnegative integer. Converting this sum into an integral over all 
positive k will remove the quantization radius L. The eigenenergies for s-states are given 
by the equation 

a cot(aR) + y£ = . 

Within these results the matrices v and u can be evaluated via Eqs. (|60| , |61| ) and all the 
theory described above can be applied in a straightforward way. 

As a realistic physical picture, we take the depth of the spherical well V = 1 GeV 2 , radius 
R = 1 fm, and a lifetime of the condensate T = 1 fm/c. Then the lowest level is at a depth 

(858 MeV) 2 which corresponds to Q = 



of £ =(858 MeV) which corresponds to Q = J\£\ — m 2 = 847 MeV. These assumptions 
are probably exaggerated as the average number of pions of one particular type in this case 
is n w 2500 . 

Figure |^ displays the behavior of the dimensionless variables J* -1 ^, I^/Cl and their 
sum /( _1 )fi + jW/fi shown in dashed, dotted and solid lines, respectively, as a function of 
the potential depth, left panel, and the potential size, right panel. The singularity of 1^ /Vt 
at threshold corresponds to Q approaching zero and does not have any physical significance 
because the approximation exp(fiT) ^> 1 does no longer hold. In general both plots are 
dominated by the behavior of f2 for the chosen parameters whereas the integrals /( -1 ) and 
1^ are weakly influenced by the form of II(x). 
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FIG. 6. The behavior of the dimensionless variables I^ 1 '^}, I^'/Q, and I^^Q + I^- 1 ' /&, is shown 
in dashed, dotted and solid lines, respectively, as a function of the potential depth W, left panel, 
and potential size R, right panel. A fixed size of 1 fm was used for the left panel and a fixed depth 
^f =1 GeV was used for the plot in the right panel. 
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FIG. 7. The average number of particles produced as a function of the depth of the potential 



field v V is shown in the left panel, the size R was fixed at 1 fm. The right panel shows the number 
of particles produced as a function of size R given a fixed depth V = 1 GeV 2 . The time length of 
the perturbation is set at T =1 fm/c. 

Fig. [7] shows the average number of particles as a function of the well size, right panel, and 
depth of the perturbation, left panel. The average number of particles grows approximately 
exponentially with the potential depth (Fig. |7|, left panel). This is simply related to the 
fact that in a deep well the ground state energy grows almost linearly with the depth having 
Q, « yV. By making the potential wider the ground state approaches the bottom of the 
well, limiting Q to a constant. This restrains the growth of particles shown on the right panel 
of Fig. [7[ One should bear in mind that the conditions considered are quite extreme and 
were used here to emphasize the character of the condensate. Practically the time scale may 
be shorter and perturbation weaker, leading to a more pre-condensate picture with much 
fewer pions. The total energy available in a heavy ion reaction may provide a guideline to 
what the perturbation II is and the realistic number of mesons produced are. 

As a final part of this analysis we consider non-condensate pion production, which may 
be the main mechanism in most practical situations. Eq. ( |74|) with the additional help of 
Eq. (|75|) results in the following form for s-wave pions from the continuum 
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k 2 ) 2 u k W K (u k + W K ) 



(77) 
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where parameter a is defined as a = \/k 2 + V, W K = J k 2 + m 2 , and Wfc = Jk 2 + m 2 are the 
total energies of the corresponding modes. As expected, the number of particles produced 
with no condensate involved is quite small. The left panel in Fig. [5] shows the average 
number of non-condensate pions as a function of the potential depth for different spatial sizes 
R =0.5, 1, 2 fm. The right panel of the same figure displays the dependence on the size R for 
various values of V . The important conclusion here is that the number of non-condensate 
pions ranges from a few up to maybe a dozen for extreme cases. This number is completely 
negligible in the presence of a strong condensate that produces hundreds of mesons. However, 
present experiments may just barely reach the point of the phase transition, and therefore 
the fraction of non-condensate pions is considerable, if not dominant. Moreover, other 
(conventional) mechanisms of pion production have to be taken into account. 
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FIG. 8. The left side shows the average number of pions of a particular type as a function of 
V , the depth of the perturbation. Curves displayed as solid, dotted and dashed lines correspond 
to the values of the radius R of 0.5, 1, and 2 fm, respectively. Plotted on the right hand side is 
the number of pions versus the radius R for values of yV of 0.5, 1 and 2 GeV as solid, dotted and 
dashed lines, respectively. 



V. SUMMARY AND CONCLUSIONS 

Our prime objective in this work was to study a mechanism of pion production in heavy 
ion collisions related to the creation of a chiral condensate and to explore how pion distri- 
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butions can signal the presence of a condensate. We studied meson production by imposing 
a pion dispersion relation specific to the medium, i.e. with a space- and time-dependent 
effective mass. 

In general, the problem of parametric excitation of the field quanta presents an interesting 
question as it is encountered in many branches of physics from condensed matter to high 
energy physics. The problem also exhibits a vast variety of solutions ranging from adiabatic 
to phase transitions and condensates. 

We have conducted an extensive study of quantum field equations of a general form 
of Eq. (Q). In our picture the parametric excitation of the field quanta is carried out by 
the externally given space- and time-dependent mass term. This term in the Lagrangian 
is quadratic in the field, and the states produced are often called "squeezed" J27| , |28| , |30| , |3l| . 
Some analogy can be drawn here from well studied linear current type terms that produce 
coherent states. However, the essential difference between coherent and squeezed states 
arises due to the fact that the latter correspond to the pairwise generation of quanta. In the 
case of pseudoscalar pions, this means that charge, isospin and parity are exactly preserved. 
We have emphasized the fact that the quantum solution can be built from the classical 
solution, and this important link was established via a canonical Bogoliubov transformation. 

With our interest lying in the direction of chiral condensate we focused our attention 
on the potential of Eq. (0) to form a condensate with fast particle production and large 
correlations when the effective mass goes through zero. Along with a general formalism 
that can be used for numerical studies we have analytically solved the problem when the 
effective mass experiences sudden abrupt changes. Consideration of this particular temporal 
perturbation allowed us to clearly separate the exponentially rising collective pion condensate 
modes for any given spatial form of the perturbation in the effective mass. This produced 
conditions where the condensate and its signatures can be seen. 

We have identified two basic channels of pion production. The first involves only a few 
discrete condensate modes and a large associated pion population. The second leads to the 
production of far fewer (non-condensate) particles with a broad phase space distribution. 
Mathematically these channels can be identified as production of mesons from bound and 
continuum states of a Schrodinger equation with a potential of the form of the perturbation 
itself, see Fig. [5| The bound states of negative energy are responsible for the characteristic 
features of the condensate. 

Numerically, the number of non-condensed pions ranges from a few up to a dozen, and 
as expected is not very sensitive to the choice of the spatial and temporal form of the 
perturbation in the effective mass. In contrast, the condensate modes have an exponential 
sensitivity to the input parameters. As the abrupt changes in the effective pion mass grow 
in strength, a critical point is reached with the appearance of the condensate mode. The 
population of this mode increases dramatically from zero to thousands with a further slight 
change in the mass parameter. Due to this hyper-sensitivity of the number of condensed 
pions to the perturbation it is practically impossible to predict the effect quantitatively 
without specifying precisely the scenario of the process. 

However, our results predict the number of non-condensate pions, thereby imposing a 
lower limit on the statistics needed to unambiguously detect the chiral condensate. Fur- 
thermore, we have shown that the condensate pions have a specific momentum distribution 
due to their common collective condensate mode. We have also shown that although the 
distribution over species starts from the famous 1/a/J form for one mode it quickly becomes 
Gaussian with the appearance of successive modes. Therefore the presence of several modes 
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complicates the detection of a chiral condensate. The number of modes present increases 
with energy. In addition, the mass parameter can be strongly perturbed in more than one re- 
gion, each increasing the number of modes. Tunneling and chaotic dynamics in the resulting 
multi-well potential might lead to another class of interesting problems. 

This work can be extended in several directions. With the formalism presented here, large 
numerical studies of the effective pionic field in a hot medium can be conducted involving 
realistic and even self-consistent forms of the perturbation given by the a field. Constraints 
on the perturbation of the mass parameter should be related more rigorously to observables. 
Further analysis of our results applied to phase transitions, zero mass particle production, 
energy transfer and many other field theory problems would definitely be fruitful. We feel 
that this work may provide a step forward in the study and classification of field theories 
with parametric excitations and possibly clarify the nature of the produced squeezed states. 
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